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We compute the electromagnetic radiative corrections to all leading annihilation processes which 
may occur in the Galactic dark matter halo, for dark matter in the framework of supersymmetric 
extensions of the Standard Model (MSSM and mSUGRA), and present the results of scans over the 
parameter space that is consistent with present observational bounds on the dark matter density 
of the Universe. Although these processes have previously been considered in some special cases 
by various authors, our new general analysis shows novel interesting results with large corrections 
that may be of importance, e.g., for searches at the soon to be launched GLAST gamma-ray space 
telescope. In particular, it is pointed out that regions of parameter space where there is a near 
degeneracy between the dark matter neutralino and the tau sleptons, radiative corrections may 
boost the gamma-ray yield by up to three or four orders of magnitude, even for neutralino masses 
considerably below the TeV scale, and will enhance the very characteristic signature of dark matter 
annihilations, namely a sharp step at the mass of the dark matter particle. Since this is a particularly 
interesting region for more constrained mSUGRA models of supersymmetry, we use an extensive 
scan over this parameter space to verify the significance of our findings. We also re-visit the direct 
annihilation of neutralinos into photons and point out that, for a considerable part of the parameter 
space, internal bremsstrahlung is more important for indirect dark matter searches than line signals. 



PACS numbers: 13.40.Ks,95.35.+d, 11.30.Pb, 98.70.Rz 



I. INTRODUCTION 

During the last few years a strong consensus has 
emerged about the existence of a sizeable dark mat- 
ter contribution to the total cosmological energy den- 
sity. The identification of experimental signatures that 
eventually may determine the nature of the cosmologi- 
cal dark matter is thus becoming ever more important. 
The present estimates [l[ give the fraction of the criti- 
cal density of cold dark matter particles as flcDAih 2 ~ 
0.105 ± 0.013, where the Hubble parameter (scaled in 
units of 100 km/s Mpc" 1 ) is h ~ 0.70 ± 0.02. Also, on 
the scales of galaxies and smaller, a number of meth- 
ods including measurements of rotation curves as well 
as gravitational lensing agree well with the predictions 
from N-body calculations of gravitational clustering in 
cold dark matter cosmologies (see e.g. @). 

The methods of detection of dark matter (for reviews, 
see Q) can be divided into accelerator production and 
detection of missing energy (especially at the LHC at 
CERN, which will start operating some time in 2008), 
direct detection (of dark matter particles impinging on a 
terrestrial detector, with recent impressive upper limits 
reported by Q), or indirect detection of particles gener- 
ated by the annihilation of dark matter particles in the 
Galactic halo or in the Sun/Earth. All these methods are 
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indeed complementary - it is probable that a signal from 
more than one type of experiment will be needed to fully 
identify the particle making up the dark matter. The 
field is just entering very interesting times, with the LHC 
soon starting and new detectors of liquid noble gases be- 
ing developed for direct detection. For indirect detection, 
the satellite PAMELA [| was launched a year ago and 
will soon reveal its first sets of data for positron and an- 
tiproton yields in the cosmic rays Q. AMANDA Q at 
the South Pole that has searched for detection of neutri- 
nos from the centre of the Earth or the Sun Q , will soon 
give way to the much larger detector IceCUBE ||, and 
for gamma-rays coming from annihilations of dark mat- 
ter particles in the halo [l(| the space satellite GLAST 
pd| . to be launched in 2008, will open up a new window 
to the high-energy universe, for energies from below a 
GeV to about 300 GcV. 

One problem with all these discovery methods is that 
the signal searched for may be quite weak, with much 
larger backgrounds in many cases. For indirect detection 
through gamma-rays, the situation may in principle be 
better, due to (i) the direct propagation from the region 
of production, without significant absorption or scatter- 
ing; (ii) the dependence of the annihilation rate on the 
square of the dark matter density which may give "hot 
spots" near density concentrations as those predicted by 
N-body simulations; (iii) possible characteristic features 
like gamma-ray lines or steps, given by the fact that no 
more energy than m x per particle can be released in the 
annihilation of two non-rclativistic dark matter particles 
(we denote the dark matter particle by \)- 

As an example, it was recently shown [l2j that in mod- 
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els of an extended Higgs sector, the line signal from the 
two-body final states 77 and Z7 could give a spectac- 
ular signature in the gamma-ray spectrum between 40 
and 80 GeV. On the other hand, in models of univer- 
sal extra dimensions (UED) [H]) or in the theoretically 
perhaps most favoured, supersymmetric, models of dark 
matter the line feature is in general not very prominent, 
except in some particular regions of the large parameter 
space. However, it was early realised that there could be 
other important spectral features [l4| , and recently it has 
been shown that internal bremsstrahlung (IB) from pro- 
duced charged particles in the annihilations could yield 
a detectable "bump" near the highest energy for heavy 
gauginos or Higgsinos annihilating into W boson pairs, 
such as expected in split supersymmetry models |15| . In 
[l6| , it was furthermore pointed out that IB often can be 
estimated by simple, universal formulas and often gives 
rise to a very prominent step in the spectrum at photon 
energies of E 1 = m x (such as in UED models [13 )■ 

Encouraged by these partial results, we have performed 
a detailed analysis of the importance of IB in the min- 
imal supersymmetric extension to the standard model 
(MSSM). We have therefore calculated the IB contribu- 
tions for all two-particle charged final states from general 
neutralino annihilations. Besides confirming the men- 
tioned partial results for the universal radiative correc- 
tions, in particular those relating to soft and collincar 
bremsstrahlung, we also point out interesting cases of 
model-dependent "virtual" bremsstrahlung (i.e. photons 
emitted from charged virtual particles), see Fig. 1. We 
confirm the suspicion expressed already in [l4| that this 
type of emission may circumvent the chiral suppression, 
i.e., the annihilation rate being proportional to to 2 for 
annihilation into a fermion pair from an S'-wave initial 
state, as is the case in lowest order for non-relativistic 
dark matter Majorana particles in the Galactic halo (see 
also [HI)- Since this enhancement mechanism is most 
prominent in cases where the neutralino is close to de- 
generate with charged sleptons, it is of special impor- 
tance in the so-called stau coannihilation region in mod- 
els of minimal supergravity (mSUGRA, as implemented 
in [3). We therefore run through an extensive scan over 
these models (based on [2(|) and find, indeed, remarkable 
cases of enhancement of the gamma-ray rate in the stau 
coannihilation region, near the maximal possible photon 
energy E 1 = m x . 

Let us stress that the radiative corrections to the main 
annihilation channels, here computed systematically for 
the first time, may turn out to be of utmost importance 
when fitting gamma-ray data, e.g. from GLAST, to su- 
persymmetric dark matter templates. Over much of the 
parameter space we have scanned, these corrections give 
a large factor of enhancement over the commonly adopted 
estimates, especially at the observationally most interest- 
ing, highest energies. More importantly, they add a fea- 
ture, the very sharp step at the dark matter mass, that 
would distinguish this signal from all other astrophysical 
background (or foreground) processes. 



II. SUPERSYMMETRIC EXTENSIONS TO THE 
STANDARD MODEL 

Although our calculated electromagnetic radiative cor- 
rections will be applicable to the annihilation of any Ma- 
jorana dark matter WIMP (weakly interacting massive 
particle), we will present results for the arguably most 
plausible dark matter candidate; the lightest neutralino 
in the MSSM, which is a linear combination of the su- 
perpartners of the gauge and Higgs fields: 

X = X°i = N U B + N 12 W 3 + N 13 H° + N 14 H° . (1) 

We perform all numerical calculations using the Dark- 
SUSY code (see [2l| for our sign conventions and other 
details). 

The parameter /1 is as usual the Higgsino mass param- 
eter, tan (3 is the ratio of vacuum expectation values of 
the two Higgs doublets, M\, M 2 and M 3 are the gaugino 
mass parameters, uia is the neutral pscudoscalar Higgs 
mass. We are also using parameters too, A t and A^, 
defined through the simplifying ansatz: Mq = My = 
M D = M E = M L = m l, A v = diag(0, 0, A t ), A D = 
diag(0, 0, Ab), Ag = 0. Here A are soft trilinear cou- 
plings and M soft sfermion masses which in general are 
3x3 matrices in generation space, but are thus simplified 
through our ansatz and encoded in too, A t andA^. We do 
not allow for CP-violating phases other than the CKM 
phase of the standard model. As a natural further sim- 
plification, the grand unification condition for the gauge 
couplings, leading to Mi = | tan 2 9 W M 2 « |M 2 is used. 
For the MSSM scans, we use FeynHiggsFast [22j for the 
Higgs boson masses and decay widths. For each model, 
we will denote with m x the mass of the lightest neutralino 

and with Z g = |An| 2 + | A^i2 1 2 the gaugino fraction. 

As a more restricted, but in some sense more natural, 
set of parameters we use those that stem from demanding 
that the electroweak symmetry be spontaneously broken 
by electroweak radiative effects such a s ap pears in mini- 
mal supergravity (mSUGRA) models 12311 . Also here we 
use the implementation in DarkSUSY [2l| . which relies on 
the public code Isajet [l|| for the solution of the renor- 
malization group equations (RGE) and for the mass spec- 
tra. In these models, parameters are given at the grand 
unification scale, and are then (using the RGEs) let to 
run down to the electroweak scale. The range of models 
with correct symmetry breaking are usually displayed in 
the TOo-mi/2 plane, where too and m.1/2 are the universal 
scalar and gaugino masses, respectively, at the grand uni- 
fication scale. The additional parameters of mSUGRA 
models are tan/3 (as for the MSSM), Aq (which is the 
common trilinear term at the grand unified scale) and 
the sign of fj, (\fj,\ is determined from the other parame- 
ters) . 

Let us briefly mention those regions in the mo-mi/2 
plane that are important from a cosmological point of 
view, as they correspond to models with a relic den- 
sity in accordance with the WMAP value: The bulk 
region which has low Too and TO1/2; the funnel region 
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FIG. 1: Types of diagrams that contribute to the first or- 
der QED corrections to WIMP annihilations into a pair of 
charged particle final states. The leading contributions to di- 
agrams (a) and (b) are universal, referred to as final state 
radiation (FSR), with a spectral distribution which only de- 
pends slightly on the final state particle spin and has been 
calculated, e.g., in fid/ . Internal bremsstrahlung from virtual 
particles (or virtual internal bremsstrahlung, VIB) as in dia- 
gram (c), on the other hand, is strongly dependent on details 
of the short- distance physics such as helicity properties of the 
initial state and masses of intermediate particles. 



niA ~ 2m x , where annihilations in the early universe 
are enhanced by the presence of the near-resonant pseu- 
doscalar Higgs boson; the hyperbolic branch or focus 
point region where mo 3> the stau coannihilation 

region where m x f» mf ; and finally the stop coannihila- 
tion region (arising when Aq 7^ 0) where m x at mj. The 
stau coannihilation region has recently been noticed to 
have favourable properties for indirect detection rates in 
antiprotons and gamma-rays 2J|. In this paper we will 
show that, in addition, there is a great enhancement of 
the high energy gamma-ray signature in this region. 



III. INTERNAL BREMSSTRAHLUNG FROM 
WIMP ANNIHILATIONS 

A. The general case 

Whenever WIMPs annihilate into pairs of charged par- 
ticles XX, this process will with a finite probability au- 
tomatically be accompanied by internal bremsstrahlung 
(IB), i.e. the emission of an additional photon in the 
final state (note that in contrast to ordinary, or exter- 
nal, bremsstrahlung no external electromagnetic field is 
required for the emission of the photon). As visualized 
in Fig. [TJ one may distinguish between photons directly 
radiated from the external legs (final state radiation, 
FSR) and photons radiated from virtual charged particles 
(which we will refer to as virtual internal bremsstrahlung, 
VIB). So, to be more specific, the IB photons will be the 
total contribution from both FSR and VIB photons. 

If the charged final states are relativistic, FSR 
diagrams are always dominated by photons emitted 
collinearly with X or X. This is a purely kinematical 
effect and related to the fact that the propagator of the 
corresponding outgoing particle, 



D(p)oc ({k+p) 2 -m 2 x ) \ 



(2) 



diverges in this situation. Here, k and p denote the mo- 
menta of the photon and the outgoing particle, respec- 
tively. The resulting photon spectrum turns out to be 



of a universal form, almost independent of the underly- 
ing particle physics model (T5 , uj\ . Defining the photon 
multiplicity as 



dN 



xx 



da 



xx^xX 7 



dx 



T xx^xx 



dx 



(3) 



where x = 2E 1 /^/~s = Ey/m x and s is the center-of-mass 
energy, it is given by [161 ]: 



dN 



xx 



a( 



dx 



<x_ T 



X 



(x) log 



»(1 - *) 



•x 



(4) 



Here, Qx and mx are the electric charge and mass of X; 
the splitting function T(x) depends only on the spin of 
the final state particles and takes the form 



^fermion(^') 

for fcrmions and 



1 + (1 - .t) 2 



boson 



(x) 



1-x 



(5) 



(6) 



for bosons. Due to the logarithmic enhancement that 
becomes apparent in Eq. FSR photons are often the 
main source for IB (note that very near the kinematical 
endpoint, x ~ 1 — m 2 x /s, it is not sufficient anymore to 
only keep leading logarithms and one can thus no longer 
expect Eq. Q to be a good approximation for the actual 
spectrum) . A prominent example where FSR in this uni- 
versal form not only dominates IB but in fact the total 
gamma-ray spectrum from WIMP annihilations, is the 
case of Kaluza-Klein dark matter [lj} . 

In general, one can single out two situations where pho- 
tons emitted from virtual charged particles may give an 
even more important contribution to the total IB spec- 
trum than FSR: i) the three-body final state XXj satis- 
fies a symmetry of the initial state that cannot be satis- 
fied by the two-body final state XX or ii) X is a boson 
and the annihilation into XX is dominated by t-channcl 
diagrams. To understand that the first case only leads to 
an enhancement of VIB, and not of FSR, we recall that 
the latter is dominated by collinear photons, i.e. the (vir- 
tual) final state particles are almost on mass-shell; the 
two- and three-body final states are thus bound to the 
same symmetry constraints. The enhancement of the an- 
nihilation rate in the second case follows from a closer in- 
spection of the t-channcl propagator. For non-rclativistic 
WIMPs, it takes the form 



D t (p) cx {(l~pf 



m 2 ~ + m x + 2m x Ex 



(7) 



where I is the momentum of one of the ingoing WIMPs 
and X denotes the particle that is exchanged in the t- 
channcl. If \ and X are almost degenerate in mass, 
one thus finds an enhancement for small Ex which - for 
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kinematical reasons - corresponds to large photon ener- 
gies Ej. Note that this enhancement is less efficient for 
fermions since the infrared external spinor leg would re- 
sult in a further suppression factor. In both cases, and in 
contrast to the situation for FSR, the resulting spectrum 
is highly model-dependent. Prominent examples where 
VIB dominates over FSR as a consequence of these two 
special situations are neutralino annihilations into lep- 
tons [HI and charged gauge bosons 15 1, respectively, to 
which we will shortly return. 



B. The neutralino case 

The relevant final states of neutralino annihilations arc 
W+W~, W ± H^, H+H- and //; analytical expressions 
for the corresponding rates can be found, e.g., in Q. 
The inclusion of an additional photon in the final state is 
straight-forward, though tedious. Here, we are interested 
in the annihilation of neutralinos today, so we work in the 
limit of vanishing neutralino velocity - which greatly sim- 
plifies the calculation and the form of the analytical ex- 
pressions. Our analytical results for arbitrary neutralino 
compositions agree, in the corresponding limits, with the 
special situations considered earlier, i.e. pure Higgsino or 
Wino annihilation into W + W~ [15[ and photino annihi- 
lation into light leptons [T3|. Let us stress that, although 
we state in the following simplified expressions for the 
photon multiplicities in some special cases, we always 
use the full expressions in our numerical calculations. 

Let us first note that for neutralino annihilations, in 
contrast to the situation for Kaluza-Klcin dark matter, 
we cannot in general expect very large contributions from 
FSR (from external legs). This is because the light- 
est charged final states, for which the logarithmic en- 
hancement shown in Eq. ([4J would be most effective, are 



fermionic - but the annihilation rate of neutralinos into 
light fermions is strongly suppressed by a factor m 2 /m x 
due to the helicity properties of a highly non-relativistic 
pair of Majorana fermions [25|. On the other hand, as we 
will demonstrate now, all the possible final states for neu- 
tralino annihilations have the potential of showing con- 
siderable VIB contributions. 

For large neutralino masses m x 3> myy and charginos 
almost degenerate with the neutralino, e.g., W + W~ and 
W final states fall into the second of the categories 
discussed in the previous subsection. If the neutralino is 
an almost pure Higgsino (or Wino) , we then find that the 
photon multiplicity for W + W~ final states can be well 
approximated by (see also Ref. [l5[ for a similar result): 



dN w+ 



w- 



4(1 



„2\2 



da; 



7r (l-x + e/2)x 
l-x + e/2 



(8) 



log 2 



- l/2 + x 



where e = my/ lm x . As demonstrated in the next section, 
the most significant IB contributions from final states 
with charged gauge bosons to the total annihilation spec- 
trum are, in fact, often of this form. The full analyt- 
ical expressions for dN^+w- /dx and dN^±^T /dx are 
rather lengthy and we will therefore not explicitly state 
them here. Let us stress again, however, that we of course 
do use them in our actual calculations. 

As mentioned before, the annihilation into light 
fermions is helicity suppressed; for large photon ener- 
gies, however, fermion final states containing an addi- 
tional photon, //7, are not subject to such a suppres- 
sion. While our full analytical expressions are again 
rather lengthy, they simplify considerably in the limit of 
nif — > 0. In this case, and assuming that both sfermions 
have the same mass, the photon multiplicity is given by 



dN 



f+f- 
dx 



— UcmV/ 



IslI 4 



64tt 2 



m 2 x (av) xx ^ ff 



x(l-*){- 



4.c 



2x 



{l + H)(l+H-2x) (1 + n-x) 2 



(1 +£)(1 2x) 

{l+H-xf 



log 



1 + f-l ■ 

1 + n - 2x . 



(9) 



where fx = m 2 - /m x = m 2 ~ /m 2 , and qrPl (ffi-ffl) is the 
coupling between neutralino, fermion and right-handed 
(left-handed) sfermion. This confirms the result found 
in [3| for photino annihilation (while [IB] states a re- 
sult that is an overall factor of 2 larger). Note the large 
enhancement factor m x /m 2 due to the lifted helicity 
suppression (from {cv) _^fj <x m 2 mx~ A ), and another 
large enhancement that appears at high photon energies 
for sfermions degenerate with the neutralino. 

Charged Higgs pairs H + H~, finally, provide yet an- 



other interesting example of the second category dis- 
cussed in the previous subsection: In the limit v —> 0, 
this final state is not allowed because of CP conservation. 
The annihilation into H + H~~/, on the other hand, is pos- 
sible. However, since charged Higgs bosons in most mod- 
els have considerably larger masses than gauge bosons the 
enhancement mechanism described in the previous sub- 
section is not as efficient as in the latter case. These final 
states are thus expected to be of less importance in our 
context. 
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In the last part of this section, let us briefly describe 
how we implemented IB from the various possible final 
states of neutralino annihilations in DarkSUSY. The total 
gamma-ray spectrum is given by 




where Bf denotes the branching ratio into the annihi- 
lation channel /. The last term in the above equation 
gives the contribution from the direct annihilation into 
photons, 77 or Zj, which result in a sharp line feature 
[271 ]. The first term encodes the contribution from sec- 
ondary photons, produced in the further decay and frag- 
mentation of the annihilation products, mainly through 
the decay of neutral pions. This "standard" part of the 
total gamma-ray yield from dark matter annihilations 
shows a feature-less spectrum with a rather soft cutoff 
at E-v 



m 



x ■ 



In DarkSUSY, these contributions are in- 
cluded by using the Monte Carlo code PYTHIA [H[ to 
simulate the decay of a hypothetical particle with mass 
2m% and user-specified branching ratios Bf. In this way, 
also FSR associated to this decay is automatically in- 
cluded (the main contribution here comes from photons 
directly radiated off the external legs, but also photons 
radiated from other particles in the decay cascade are 
taken into account). On the other hand, IB from the 
decay of such a hypothetical particle cannot in general 
be expected to show the same characteristics as IB from 
the actual annihilation of two neutralinos. In particular, 
and as discussed in length at the beginning of this Sec- 
tion, we expect important VIB contributions in the latter 
case - while in the first case there are simply no virtual 
particles that could radiate photons. We therefore calcu- 
late analytically the IB associated to the decay (i.e. FSR 
from the final legs) and subtract it from dNj' scc /dx as 

obtained with PYTHIA; for dN]' m /dx, we then take the 
full IB contribution from the actual annihilation process 
as described before. Hence, this procedure leaves us with 
corrected PYTHIA results without FSR on the external 
legs and our analytical calculation of IB (including FSR 
and VIB) that wc add to this. 1 

Let us conclude this section by showing in Fig. [3] four 
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1 We would like to stress that this prescription is fully consistent 
since both the original and the corrected IB versions are gauge- 
invariant separately. Strictly speaking, however, we have only 
corrected for photons originating directly from the external states 
and not for those radiated from particles that appear later in the 
decay cascade. On the other hand, one would of course expect 
that modifying the energy distribution of the charged particles 
corresponding to these external legs also affects the further de- 
cay cascade. Note, however, that the resulting change in the 
photon spectrum is a second order effect; more important, for 
kincmatical reasons it does not affect photons at energies close 
to m x - which, as we shall see, are the most relevant. Finally, we 
observe that our subtraction procedure has only a minor effect 
on the photon spectrum obtained by PYTHIA and no practical 
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FIG. 2: From top to bottom, the gamma-ray spectra for the 
benchmark models defined in Tab. [T] is shown. The contribu- 
tions from IB and secondary photons is indicated separately 
(in these figures, the line signal is not included). 
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BM4 


2.2-10 4 


7792 


24.1 


17.7 


+1 


1926 


1.2-10~ 4 


0.11 


X^(1996) 


0.012 


10.8 


2.1 



TABLE I: Benchmark models that represent typical regions in the supersymmetric parameter space where IB becomes impor- 
tant. The "t-channel" entry indicates the main contributing t-channel diagram, with the corresponding sparticle and its mass 

(all masses are given in GeV). S = N-, 10 --ii° ( ^l:i !i -i ( iooGeV ) ^ s ^ Ile resca l e< i flux from IB alone and the last two columns give 
the ratio of the integrated flux, all above 0.6 mx, between the new IB contribution and secondary photons as well as the line 
signals. The main difference between BM2 and BM3 is that in the former neutralinos mainly annihilate into r leptons, while 
in the latter mainly into t quarks. Also, for BM2 only the r final states give an important contribution, while in the second 
case, even the other leptonic final states contribute considerably (due to near-degenerate slepton masses) . For the BM4 model, 
the IB contribution from the W + W~ state dominates. 



typical examples of mSUGRA models with particularly 
pronounced IB features. From top to bottom, they show 
situations in which IB from tt, t + t~ , all lepton and 
W + W~ final states, respectively, dominates the total 
gamma-ray spectrum from DM annihilations. 



IV. NUMERICAL RESULTS 

We are now ready to apply the results in the previous 
section to a set of supersymmetric MSSM and mSUGRA 
models. For MSSM, we use a set of scans of the 7 param- 
eters mentioned in Sec.|TTJ These scans are fairly general 
and partly use the method of Markov Chain Monte Carlo 
models to focus on models that have the relic density 
in the WMAP preferred range. In total we have about 
200 000 MSSM models that pass the WMAP relic den- 
sity constraint and all accelerator constraints (checking 
experimental bounds on various masses and branching 
ratios, such as that of b — > 37). For mSUGRA, we use 
the very large set of scans performed in Ref. [2(| (also 
made with a Markov Chain Monte Carlo). This data 
set contains about 550 000 models that pass all the con- 
straints. 

Since the IB photons are most distinguished at higher 
energies (as seen in Fig. we will integrate the flux 
above 0.6m x and compare the flux in this energy range 
from the flux obtained from secondary photons (arising 
mainly from 7r° decays in quark jets). We will also com- 
pare it with the monochromatic gamma ray lines [27j that 
appear at loop-level. 

In Fig. [3] we show our combined result for both the 
MSSM and mSUGRA models. In the left panel we show 
how IB compares with the regular secondary photons 
from quark jets. We plot this in the Z g j (1 — Z g ) vs m x 



consequences for the quantities we calculate in the next Section 
in order to quantify the importance of IB for neutralino annihi- 
lations. 



plane where gauginos (binos and winos) are at the top, 
Higgsinos are at the bottom and mixed neutralinos arc 
in the middle. In the right panel, we show how IB pho- 
tons compare with the monochromatic lines. It can be 
clearly seen that in large parts of this plane, IB photons 
are the dominating component, outnumbering both sec- 
ondary and monochromatic photons. Hence, this effect 
is very important to include when searching for gamma 
ray signatures from dark matter. 

In Fig. [3] we looked at the ratio of IB to other gamma 
ray contributions, let's now turn to the absolute fluxes. 
We can write the flux observed at earth as 



%{E,Aft, V) 



(11) 



dS 



9.35 • 10~ 14 — x ( J (ip) ) (Ail) era- 2 s" 1 sr" 1 . 
dE 

where we have put all the particle physics into the quan- 
tity S and the astrophysics into J. S is given by 



dE J 



IB 7 



100 GeV 



E 



V(Tf 



10 -29 cm 3 g-1 



dN]< IB 
dE 



(12) 

where the sum is over all final states /. The astrophysical 
part (that depends on the chosen halo profile) is given by 



;j(V0>(Afi) = 
1 1 





/ dL { 


/An 


J line \ 



of sight 



p{L, V) 
0.3GcV/cm 2 



(13) 
2 



For example, for a standard NFW [29|] halo profile (J) ~ 
21 averaged over Aft = 1 sr towards the galactic centre. 
This value can be orders of magnitude larger though both 
for steeper profiles and/or for more concentrated obser- 
vations towards the galactic centre (i.e. smaller AJ7 - as, 
e.g., for Air Cherenkov Telescopes). The effect of specific 
halo profiles and a more detailed analysis of the absolute 
fluxes with respect to the astrophysical background will 
be left for future work [30j , We will here focus on the 
particle physics factor S. 
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FIG. 3: Integrated internal bremsstrahlung flux from supersymmetric dark matter, above 0.6 m x , as compared to the "standard" 
continuum flux produced by secondary photons (left) and the flux from both line signals (right). As for the following figures (|4] 
and[5}, two symbols at the same location always indicate the whole interval between the values corresponding to these symbols. 
Every model considered here features a relic density as determined by WMAP and satisfies all current experimental bounds. 
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FIG. 4: The observationally relevant quantity S 
and right panel). See text for more details. 



N. 



^10- 29 cm 3 s- 1 VlOOGcV/ 



for IB (left panel) and the line signals (middle 



In Fig.|l]we show the quantity £, which is dS/dE inte- 
grated above 0.6 m x . In the left panel, we show the yields 
S for the IB contribution, in the middle for monochro- 
matic 77 and on the right for Zj. In the regions where 
the IB contribution was the largest in Fig. O we typi- 
cally have lower absolute yields. However, there are very 
pronounced regions, especially at small and intermediate 
masses, where the IB yields are very high even in ab- 
solute terms. We also note that, for ncutralino masses 
in the TeV range, we expect a sizeable increase of the 
annihilation rate due to non-perturbativc effects related 
to long-distance forces between the annihilating particles 
[3lT | . These effects have not been taken into account here 
and would result in a considerable enhancement (by a 



similar factor) of the quantity S for both line signals and 
IB. 

In Fig. \E\ we focus on the mSUGRA case and show the 
contribution relative to the secondary yield of gamma 
rays for various final states separately. In the left panel, 
we show the IB yield from the W + W~ channel, in the 
middle from the t + t~ channel and in the right from the 
ti channel. Large IB contributions for the W + W~ chan- 
nel occur when a chargino is almost degenerate with the 
neutralino, as is the case for the focus point region. Note 
that due to the grand unification condition, M\ f=a ijW 2 , 
a large gaugino fraction Z g always means that the neu- 
tralino is a Bino, with vanishing annihilation rates to 
W + W~ or W + W~j final states. The large yields from 




FIG. 5: As in the left panel of Fig. [3] but now for the individual contributions from various final states of neutralino annihilations 
in mSUGRA models. IB from light leptons covers a very similar region of the plotted parameter space as that from r leptons. 
All (other) final states not shown here give always IB fluxes less than 10% of the flux from secondary photons. 



the £ + £~ and tt channels, on the other hand, occur when 
there is a strong degeneracy with the lightest £ and i re- 
spectively. These latter cases occur in the phenomcnolog- 
ically important f and t coannihilation regions: in these 
regimes, coannihilations with f and t, respectively are 
needed to push the relic density down into the WMAP 
preferred region. Hence, we have a strong mass degener- 
acy between \ an d f/t which forces the IB contribution 
to the gamma yields to be strong. 

As for the other possible final states, we note that the 
corresponding IB contributions never exceed 10% of the 
secondary photon flux; these channels are subdominant 
also for the MSSM models contained in our scan. In 
fact, from our discussion in the previous section, this is 
somewhat expected: Charged Higgs bosons, for example, 
are always heavier than charged gauged bosons, so multi- 
TeV neutralino masses would be needed for sufficiently 
large annihilation rates into W^H^j or H + H~j (recall- 
ing that the annihilation rate in these cases is enhanced 
for relativistic final states). IB from light quarks is sup- 
pressed by the mass difference between the neutralino 
and the corresponding squark (as compared to the small 
mass difference that can be achieved in the stop coan- 
nihilation region); down- type quarks, finally, receive a 
further suppression due to their smaller electric charge. 

The main results of our paper may be more easily 
grasped by looking at the effect of IB on a small number 
of benchmarks models. Of course, for the mSUGRA case, 
it is known that the exact location in parameter space of 
such benchmarks depends very sensitively on details of 
the calculation (see e.g. [Hf). We therefore define our 
own set in Table HI which is very similar to that used 
by [2(i| except that we also include one point in the focus 
point region (BM4). This set of benchmark models is cal- 
culated with ISAJET 7.69 together with DarkSUSY 
(see [2(| for details). Point BM1 is a model where Aq has 
been chosen large and negative to make the stop almost 
degenerate with the neutralino. BM2 is a model where 



the stau is almost degenerate with the neutralino and 
in BM3 also the selectron and the smuon are degenerate 
with the neutralino. BM4, finally, is in the focus point re- 
gion, i.e. where the lightest chargino is almost degenerate 
with the lightest neutralino. The main IB characteristics 
of these benchmark models are summarized in Table HI 



V. SUMMARY AND CONCLUSIONS 

As can be seen already from our benchmark points 
in Table U and in more detail from the scatter plots in 
Fig. [31 the internal bremsstrahlung effects computed in 
this work can be very significant, changing sometimes 
by more than an order of magnitude the lowest-order 
prediction for the high-energy gamma-ray signal from 
neutralino dark matter annihilation. Although some of 
these enhancements have been found before [3, [TH, [l(| , 
this is the first time the first-order radiative corrections 
have been computed systematically, for all relevant final 
states in supersymmetric dark matter models. The re- 
sulting enhancements of the expected fluxes are surpris- 
ingly large over significant regions in the parameter space 
of the MSSM, including the more constrained mSUGRA 
models. Despite the fact that some large corrections ap- 
ply to absolute rates that are too small to be of practical 
interest, Fig. 4 shows that the quantity S, which is di- 
rectly proportional to the expected signal in gamma-ray 
detection experiments, also is significant for the internal 
bremsstrahlung contribution in large regions of param- 
eter space. For m x < 300 GcV, for example, values of 
Sib greater than 0.1 are generic, and for masses below 
100 GeV, values of 1 or higher are common, which in 
very many cases is higher than the corresponding values 
for the line signals 77 and Zj. One should also bear 
in mind that the sensitivity of Air Cherenkov Telescopes 
increases significantly with energy; detectional prospects 
for a m x ~ 1 TcV neutralino with S ~ 0.01, e.g., cor- 
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respond very roughly to those for a m x ~ 100 GeV 
neutralino with S ~ 0.5 (see, e.g., [11]). In this light, 
the situation becomes very interesting even for TeV scale 
Higgsinos, where IB generically contributes more than 10 
times as much as secondary photons. 

We note that (as anticipated in [3]) helicity suppres- 
sion and also CP selection rules of certain final states 
may be circumvented by emitting a photon; this is for 
example the origin of the very substantial enhancements 
of the signal obtained in the stau annihilation region 
in mSUGRA models. In this situation, the probability 
of emitting gamma rays vanishes at zero photon energy 
but increases rapidly at high energy (see Fig. (2}d and 
Gfc), which gives a photon "bump" at E 1 sw m x . We 
also note that in less constrained versions of the MSSM 
than considered here, we expect even more situations 
where large enhancements of the annihilation signal due 
to internal bremsstrahlung can be found. An example is 
heavy Wino dark matter [HI, Hi]], which becomes possible 
when relaxing the condition Mi as \M<z (as realized, e.g., 
in anomaly mediated supersymmetry breaking scenarios 

Hi). 

Of course, the line signals, in particular 77, have the 
virtue of being at the highest possible energy, so in order 
to make a more accurate comparison between these and 
the IB signal computed here, one would have to model 
also the expected spectral shape of possible astrophys- 
ical gamma-ray backgrounds and the energy resolution 
of the detector. This is left for future work [3(| . We 
note, however, that in general also the new contributions 



have a characteristic signature, (c/. Fig. 2) which can 
hardly be mimicked by any known astrophysical gamma- 
ray source. In fact, in some cases these spectra could 
even be used by future experiments to distinguish be- 
tween different dark matter candidates (note that, e.g., 
the distinction between Kaluza-Klein dark matter and a 
neutralino in the focus point region like our BM4 point 
would be possible already with the energy resolution of 
present Air Cherenkov Telescopes [Hj]). 

To conclude, we have shown that the commonly ne- 
glected first-order radiative corrections to neutralino 
dark matter annihilation should definitely be taken into 
account when predicting rates for gamma-ray telescopes. 
In particular, the soon to be launched GLAST space tele- 
scope pH will have an enhanced possibility over what 
has previously been assumed to detect radiation from 
supersymmctric dark matter annihilation. The routines 
needed to compute these new processes will be included 
in the next release of the DarkSUSY package pll lo7l]. 
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